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Summary 

This report describes analytical methods for calculating 
stresses and damage caused by degradation of the matrix 
constituent in polymer matrix composite materials. Laminate 
geometry, material properties, and matrix degradation states 
are specified as functions of position and time. Matrix shrink- 
age and property changes are modeled as functions of the 
degradation states. The model is incorporated into an existing 
composite mechanics computer code. Stresses, strains, and 
deformations at the laminate, ply, and micro levels are calcu- 
lated, and from these calculations it is determined if there is 
failure of any kind. The rationale for the model (based on 
published experimental work) is presented, its integration into 
the laminate analysis code is outlined, and example results are 
given, with comparisons to existing material and structural 
data. The mechanisms behind the changes in properties and in 
surface cracking during long-term aging of polyimide matrix 
composites are clarified. High-temperature-material test meth- 
ods are also evaluated. 


1.0 Introduction 

Polymer matrix composite materials are increasingly being 
used for moderately high-temperature (200 to 350 °C) appli- 
cations. These applications do not have extreme temperature 
requirements, but still involve service temperatures well above 
the capabilities of typical epoxy matrices in current aerospace 
use. Advanced polyimide matrices have been formulated for 
service at these temperatures, and composites employing these 
matrices are already in use in engine cowls, ducts, thrust 
reversers, and such. These composites have been identified as 
an enabling technology for advancements in turbine engines, 
hot aircraft structures, and high-speed aircraft structures. How- 
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ever, a critical problem preventing more widespread use of 
these materials is the aging and degradation of the matrix. 
Continuous exposure to high temperatures has been shown to 
cause mass loss, degradation of properties, shrinkage, and 
cracking. Our inability to quantify matrix degradation and its 
effect on composite properties and behavior is a limiting 
factor in the application of composite materials to these 
problems. 

A better understanding of the degradation of high- 
temperature polymeric composite materials is essential. 
Obviously, degradation is a complex and coupled process: 
heat, moisture, and oxygen are transported into the material; 
the matrix material changes chemically; the changes in the 
matrix affect the behavior of the composite material at the 
ply and laminate level; and this behavior sometimes leads to 
cracking or other failures. Failures such as matrix cracks may 
in turn lead to greatly enhanced oxygen ingress, greater 
material changes, and more cracking — clearly a coupled, and 
potentially unstable, situation. We need to be able to analyti- 
cally model this problem. Testing is also required, but it 
cannot be pursued blindly. There must be an understanding 
of what a given test is measuring. Even without a full set of 
material property data, analytical tools can aid in the identifi- 
cation of critical phenomena, the design of tests, and the 
interpretation of test results. Coupled with experimentally 
collected material property data, analytical tools will ulti- 
mately provide design capability for durable high-temperature 
composite structures. 

The work presented here concentrates on a current weak 
point in analytical capability — linking matrix material degra- 
dation and shrinkage to the ply and laminate behavior that 
results. The objective of this report is to describe a new set of 
analytical methods that can be used to predict the effects of 
matrix degradation and shrinkage on the performance of com- 
posite laminates. 
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2.0 Background 

Extensive testing has been carried out to characterize the 
performance of polymer matrix composites at high tempera- 
tures. This report will concentrate on the most commonly 
used material systems, the NASA-Lewis-developed polymer 
monomeric reactant (PMR) polyimide materials. Other mate- 
rials, notably other polyimide thermosets and thermoplastic 
materials such as PEEK, have also been included in aging 
studies. Many of the phenomena observed are similar to 
those reported for PMR materials. Recent developments in 
the PMR material system are discussed in reference 1. Most 
testing has concentrated on isothermal aging, that is, the 
effects of continuous soaking at a single high temperature 
(refs. 2 to 10). 

Exposing PMR matrix materials to high temperatures for 
moderate lengths of time (hundreds of hours) first leads to 
increases in the glass transition temperature T . This change is 
thought to be due to additional chemical cross-linking in the 
matrix. In the presence of air, oxidation may speed this pro- 
cess (ref. 2). Other changes include increases in the material’s 
stiffness, bending strength, and shear strength, when these 
properties are measured at the exposure temperature. When 
they are measured at room temperature, these properties, as 
well as interlaminar fracture toughness, decrease (refs. 3 
and 4). The viscoelastic behavior is also affected since com- 
pliance decreases with exposure time (ref. 5). 

At the upper end of the temperature range or at longer 
exposure times, severe mass loss is also noted. Although 
fibers can lose mass, the mass loss in the composite is pre- 
dominantly matrix mass loss (ref. 6) that is caused by oxygen 
cutting up the polymer and ultimately destroying it. In non- 
reactive atmospheres, mass loss is reduced but not entirely 
eliminated (ref. 2). Mass loss is associated with notable shrink- 
age of the material and a decline in the values of the material’ s 
properties. The degradation begins at exposed surfaces, 
where a distinct degraded layer is observed in both neat resin 
specimens and composites (ref. 7). As aging continues, the 
degraded layer grows. This process is thought to be controlled 
by the rate of oxygen diffusion into the material (R.C. Hipp, 
D.M. Hannon, and P.S. McClellan, 1994, MacDonnell 
Douglas Co., unpublished data). 

Eventually, voids and cracks form in the degraded region. 
Cracks allow oxygen to penetrate deeper into the material, 
thereby accelerating the degradation process (ref. 7). Bowles 
(K.J. Bowles, NASA Lewis Research Center, unpublished 
data) suggested that in extreme cases an interconnected net- 
work of microcracks eventually forms, which allows 
extensive oxidation throughout the laminate and causes the 
composite to become completely ineffective. 

Attempts have been made to analyze and model various 
aspects of this problem. Mass-loss rates have been empiri- 
cally fitted to Arrhenius rate curves (M.A. Grayson and C.G. 
Fry, 1994, MacDonnell Douglas Co., unpublished data). More 


sophisticated models have combined modeling of oxygen dif- 
fusion into the material with the chemical reaction rate equa- 
tions to predict mass loss and the growth of degraded layers 
(R.C. Hipp, D.M. Harmon, and P.S. McClellan, 1994, 
MacDonnell Douglas Co., unpublished data). Cracking of 
these degraded layers has been predicted by a finite element 
model (ref. 4) that treated the degraded layers as layers with 
different material properties. In addition, changes in viscoelas- 
tic properties have been incorporated into standard viscoelas- 
tic analysis techniques (ref. 5). 

More general analysis techniques that can be applied to 
some aspects of this problem are also available. Although 
there are many such techniques, here only those used in this 
work are presented. Given the matrix and fiber properties 
(ref. 1 1), existing micromechanical models can predict com- 
posite properties and behavior. The effects of temperature and 
moisture on composite properties can also be predicted 
(refs. 12 and 13). Perhaps existing semi-empirical analysis 
methods for composite durability (ref. 14) could be applied to 
this problem. Success, however, is unlikely, since these meth- 
ods do not account for the unique physics observed in the 
high-temperature degradation process. 

No current analysis technique can predict the properties of, 
and stresses in, a general composite laminate undergoing 
high-temperature degradation. Changes in matrix properties, 
mass loss, shrinkage and void formation — and the uneven 
distribution of these degradation effects — all must be accounted 
for in such an analysis. 

3.0 Statement of Problem 

Given the geometry of a laminate, the applied mechanical 
loads, the constituent material properties, and the tempera- 
ture, moisture, and matrix degradation states for each ply in 
the laminate, we must find the stresses, strains, and deforma- 
tions at the laminate, ply, and microstructural levels, and 
determine if any failures occur. Once this capability is achieved 
and added to an existing composite analysis package, we must 
verify it by comparing the results to existing experimental 
data. 


4.0 Approach 

The problem is approached herein in three steps: 

(1) Based on existing experimental evidence, an empirical 
model is developed in which shrinkage and changes in matrix 
properties are functions of several different degradation 
mechanisms. 

(2) From fundamental micromechanics principles, the ef- 
fects that these changes in the matrix will have on unidirec- 
tional plies is predicted. 
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(3) By using modified classical -laminated-plate relations, 
the effect of ply changes on laminate behavior is predicted. 

These three models are then used to calculate stresses due 
to degradation and mechanical and hygrothermal loading at 
the laminate, ply, and micromechanical levels. Composite 
failure criteria are used to determine if failures occur. 


The metrics are defined such that they have the value 0 in 
virgin material and 1 in fully degraded material. They are 
defined as follows: 



5.0 Matrix Degradation Modeling 


m o 


( 2 ) 


5.1 Empirical Basis 

A quantitative and computationally useful model of the 
matrix shrinkage and the changes in material properties caused 
by degradation was constructed. It was based on information 
obtained in the current literature and in extensive discussions 
with the materials communities at NASA’s Lewis and Lan- 
gley Research Centers. The matrix was modeled as under- 
going two time-dependent processes, and the fibers were 
assumed to be stable. All symbols used herein are defined in 
appendix A. 

The first process is a matrix mechanism that is partly an 
extension of the cure process, in which the matrix changes by 
chemical and physical aging. This mechanism, which de- 
pends on temperature and time, is associated with small 
amounts of matrix shrinkage and mass loss (less than 1 per- 
cent) and with density and material property changes. In 
thermoset materials, this mechanism is dominated by changes 
in chemical cross-linking, so the mechanism will be referred 
to here as the cross-linking mechanism. 

The second process is a degradation mechanism due to 
oxidation and/or another mechanism that breaks down the 
polymer structure of the matrix. Temperature and oxygen 
availability govern this mechanism. In thick parts of the ma- 
trix, this mechanism appears to be diffusion limited and, thus, 
occurs mostly on the surface, causing material property changes 
(mostly degradations) and mass loss. Mass can be lost by 
erosion of the surface (making the part smaller), formation of 
voids, reduction in density of the matrix, or shrinkage of the 
matrix. In practice, density changes are not observed, whereas 
some void formation and major shrinkage are observed 
(ref. 5). This second mechanism will be referred to as the 
mass-loss mechanism. 

5.2 Degradation Metrics 

The degradation state at a point will be represented by two 
dimensionless metrics: a , the cross-linking mechanism; and b, 
the mass-loss mechanism. The first, a, is a measure of chemi- 
cal cross-linking and is taken to be proportional to the in- 
crease in the measured glass transition temperature of the 
matrix material. The second, b , is a measure of the breakdown 
of the polymer and is taken to be proportional to mass loss. 


where T is the measured glass transition temperature, m is 
the mass in a small volume of matrix surrounding the point at 
which the degradation is specified, and the subscripts o and / 
represent the initial and final states. Note that a is positive for 
an increase in T , whereas b is positive for a decrease in mass; 
the final state in the mass-loss case occurs when the mass is 
zero. These metrics specify the degradation state at a single 
point and at a given time. However, because the degradation 
state within a specimen is time dependent and nonuniform, 
these metrics must be specified as functions of time and 
position. 

53 Matrix Property Changes 

The effects of aging on the matrix’s room temperature 
properties are modeled empirically as functions of the two 
degradation metrics. The intent of these models is not to 
predict property changes, but to reduce the existing data to a 
computationally convenient form. The engineering stiffness 
properties E m and G m of the matrix are altered by aging as 
follows: 


P RT ~ P RTo 


1 + A E a (a)" 


■]-[l + A E b (b) n 


(3) 


where P RT represents the modified property used in all subse- 
quent calculations, P RTo represents the initial value of that 
property, and A E q and A E b represent changes in E m associ- 
ated with the first and second mechanisms respectively. The 
Poisson’ s ratio is assumed to be unaffected, and A E and n £ are 
determined by fitting equation (3) to available data for each 
reaction. The additional subscripts a and b refer to the reac- 
tion type defined above. 

Similarly, the room temperature matrix strengths in ten- 
sion, compression, and shear (S mT S mC , S mS ) are modified as 
follows: 


Prt = P RJo ■ l 1 + ] • [l + A S b V) n » ] (4) 

And the matrix thermal expansion coefficient a is modified 
by 
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(5) 


m o> 4bi Po 


P RT ~ P RTo 


[ 


l + Aa a (a) 


1 + A <x b (b) 


ab 


where again AS , A a, n s , and n fl are fit to available data. 

The effects of temperature on matrix properties are mod- 
eled by using an existing empirical relation (ref. 12). Strength 
and stiffness properties at temperature T are 


P(T)=P, 


T -T 
8 


\l/2 


RT 


IV 


l RT 


(6) 


where T rt is room temperature. Thermal expansion coeffi- 
cients are also modified by equation (6), but with the 1/2 
power replaced by -1/2 (from ref. 12). Degraded properties at 
temperature T are found by first applying equations (3) to (5) 
and then applying equation (6). 

5.4 Shrinkage and Void Formation 

Consider an infinitesimal cube of matrix material (fig. 1) 
with initial side length / and porosity <p Q . The cube is unre- 
strained and effectively homogeneous, with the porosity uni- 
formly distributed. The solid material (not including voids) 
within this cube has volume V f density p o , and mass m o . The 
volume of matrix material in the cube is then 

< 7 > 

and its mass is 

In response to degradation, the volume can (1) lose mass, 
which gives a new mass m; (2) shrink, which results in a new 
side dimension /; (3) develop voids (increase porosity), which 
yields a new porosity <p\ and (4) change density, which results 
in a new density p. The new volume and mass are thus 

v = {l -$)/ 3 (9) 



Figure 1 . — Cube of matrix material undergoing shrinkage, 
mass loss, and changes in porosity and density. 



Combining equations (8) and (10) 


m p (l~4>)/ 3 

m o Po (I - 0oK 3 


( 12 ) 


(13) 


m = p(l - <())/ 3 


(10) and substituting equation (12) into equation (13) gives 


A matrix shrinkage strain £* is defined as follows for a small 
strain: 


m 

m o 


p (1-0) 

Poi'-t) 


(1 



(14) 


* 




(ID 


Assuming small strains and discarding higher order terms in 
strain reduces equation (14) to 
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m p (l~0) 

m o Po (l“0 o ) 



(15) 


As discussed previously, degradation is assumed to consist of 
two primary mechanisms, both of which can result in shrink- 
age of the matrix. The shrinkage due to cross-linking is desig- 
nated £*, and the shrinkage due to mass loss is eg. Thus 


The fraction of mass loss due to changes in porosity is 
defined as an empirical parameter D. It can be determined by 
fitting to data from appropriate experiments. In terms of 
equation (21) 


D = 


0-0 o 

-3el+<p-<p 0 


( 22 ) 


£ 


m 


* 



+ 4 


5.5 Shrinkage Due to Cross-Linking 


The rest of the mass loss must be accounted for by shrink- 
age. Combining equations (21) and (22), and rearranging to 
express shrinkage and porosity changes, yields both 


The cross-linking mechanism is assumed to cause no mass 
loss or change in porosity. The matrix shrinkage due to cross- 
linking, if any, can be measured experimentally and fit to the 
following empirical relation: 

£* a = -Ca n (17) 


* 0 ~D)b 


(23) 


0 = 0 O + Db 


(24) 


where C and n are fit to data. In the absence of mass loss, jj Q Effects at the Ply Level 
porosity growth, and mass-loss-related shrinkage, equa- 
tion (15) requires that 6.1 Volume Fractions 


P = 


Po 

M £ o) 


(18) 


5.6 Shrinkage Due to Mass Loss 

The mass-loss mechanism is assumed to cause matrix shrink- 
age and increased porosity. Limited experimental evidence 
(ref. 15) indicates the remaining material has constant den- 
sity. If we assume constant density, and no shrinkage due to 
cross-linking, equation (15) can be expanded to 


m (l + ^-0-30£fc) 

m o (l ~t) 


(19) 


Consider a composite with an Initial matrix volume frac- 
tion k , an initial void volume fraction k , and an unchang- 
ing fiber volume fraction k ~ The effects of degradation on the 
volume fractions are calculated here. By definition. 


k f + k mn + k = 1 

j mo vo 


(25) 


If we assume that the voids are only in the matrix, the 
matrix porosity used in the previous section can be defined as 


Po = 


k -I- 
vo mo 


From equations (25) and (26) 


(26) 


Combining equations (2) and (19) and rearranging yields 


Ko 

=('-*/) 

k 

(27) 

Mo) 

and 

k mo= 1 


-k) 

(28) 

Expanding equation (20) as a Taylor series about $ = 0 and 
eliminating higher order terms in 0, and Ej* (which are 

assumed to be small) yields 

Similarly, 

and 

V 

-(■-*/) 

10 

(29) 

b = -3£* h +<i)-<p 0 (21) 


*m=l 

('-*> 

- 0 ) 

(30) 
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And combining equations (24), (28), and (30) yields the modi- 
fied matrix volume fraction k: 

m 

k m =l‘mo-('- k f ) Db ( 32 ) 

6.2 Ply Properties 

Ply properties are calculated from the fiber properties and 
the matrix properties as modified by aging (see eqs. (3) to 
(6)). The micromechanical relations of reference 12 are used. 
In the notation of reference 12, the subscripts /, m, and / 
denote a ply property, a matrix property, and a fiber property, 
respectively. Numerical subscripts refer to directions as shown 
in figure 2. The relations are repeated here for convenience: 



Figure 2. — Ply coordinate system. 



(40) 


63 Ply Shrinkage 

The matrix shrinkage associated with degradation causes 
deformations at the ply level. The fibers are assumed to be 
dimensionally stable, and hence restrain the shrinking of the 
matrix. This behavior is analogous to the thermal expansion 
of a ply consisting of a fiber and a matrix with different 
thermal expansion coefficients. It results in a ply deformation 
that is intermediate between the deformation of unrestrained 
matrix material and the zero deformation of bare fiber, as 
shown in figure 3. Micromechanical relations must be used to 
determine the net ply-level deformations and the residual 
stresses in the fibers and matrix. 

The micromechanical relations are analogous to those used 
to determine the ply thermal expansion coefficients. They are 
derived by modifying the micromechanics of reference 12. 
Matrix shrinkage coefficients and % m are defined for the 


Shrinking matrix 



Shrinking ply 


Stable fibers 

Figure 3. — Micromechanics of ply shrinkage. 
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cross-linking and mass-loss mechanisms, respectively, such 
that 


hence 


* r 

£ a=£m a 

4 = Xm b 


£ = — 
Sm a 


x = — 

i 


Combining equations (17) and (43) yields 
£ = -Ca (n-1) 


(41) 

(42) 

(43) 

(44) 

(45) 


For n = 1, which implies a linear relationship between 
shrinkage and this degradation mechanism, 


£ m =-c 

Combining equations (23) and (44) gives 

_ (l-p) 


(46) 


(47) 


The fibers are assumed to be stable, so all fiber shrinkage 
coefficients are zero; thus 

£/l 1 = £/22 = £/22 = ^/1 1 = X f22 = X /33 = 0 

where the subscripts are defined in the same fashion as for the 
thermal expansion coefficients (eqs. (39) and (40)). 

The micromechanics of this problem are identical to those 
of the thermal expansion of a composite consisting of fiber 
and matrix with different thermal expansion coefficients. Re- 
lations for the shrinkage coefficients of the plies can thus be 
developed in the same way as those for ply thermal expansion 
coefficients. This subject is well covered in the literature, so 
the development will not be repeated here. Instead, by direct 
analogy to equations (39) and (40), the ply degradation coeffi- 
cients are written as 


k E 

r _ m m r 
VI 1 “ F 

VI 1 


r _ r nr~\( * , kf v m^n i 

£122 ~ £m ~ V */ )[ 1 + E 


i ) 


(48) 


(49) 


/c £ 

v _ m m y 
AlW r X m 
VI 1 


nr \ . , k f v m E n 1 
X(22 - X(33 - 1 1 V V ) + F 

\ v n E m 


( 50 ) 


\x m (51) 


For n = 1 , equation (46) can be combined with equations (48) 
and (49) to yield 


r k m E m r 

9m- E c 

c m 


(52) 


and 


£(22 ~ £(33 - -yftf + 


\ 


^/ V m £ /ll 

E n\ ) 


(53) 


Similarly, equation (47) can be combined with equations (50) 
and (51) to yield 


Xn\ 


k f E m 


J (\\ 


0--P) 


X(22 ~ Xm 




(1 -D) 


(54) 


(55) 


-m j J 

The in-plane ply shrinkage strains due to degradation, «?u 
and ^ 22 ’ can now ca l cu l ate d as 


€ m -$in a + %m b 

e t22 = ^t22 a + %l22 b 


(56) 

(57) 


6.4 Microstresses 


The mismatch between the shrinkage of the matrix and the 
fibers creates stresses at the micro scale, even in plies that are 
stress-free at the ply level. By modifying the relations of 
reference 13, the microstresses due solely to the degradation 
effects can be calculated as follows; 


( £ m £ m) E m 

(58) 

- Ej- n e ni 

(59) 
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°m22A ~ i e ?22 £ m ) E m 

(60) 

1 - 

°m22B ~ G m22A 

(61) 

°f22 = °m22B 

(62) 

°m33A = a m22A 

(63) 

a m33B ~ a m22B 

(64) 

CT /33 1 a f 22 

(65) 


where / indicates stress in the fibers; m , stress in the matrix; 
and A and B , different regions within the matrix (see ref. 1 3 
for a complete description of this scheme). 

7.0 Effects at the Laminate Level 

7.1 Laminate Theory 


where T. is the transform matrix defined in appendix B, T 
indicates matrix transpose, and -T indicates inverse trans- 
pose. The laminate response, effective forces, and moments 
are cal-culated and applied in calculations analogous to those 
used in calculating thermal responses. The effective forces 
and moments on a laminate with N plies, N* and M* respec- 
tively, are 

n*=£q£( 0 (z|-*?) < 68 > 


= < 69 > 

1=1 L J 

Here z* and represent the distance from the top and bottom, 
respectively, of ply i to the laminate centerline, mdQ i is the 
rotated reduced stiffness matrix of the ply, the calculation of 
which is described in appendix B. The total laminate strains 
and curvatures, including mechanical, thermal, and hygral 
effects are 

e o = a(N + N r + N" + N* ) + b(M + M 7 ” + M h + M* j (70) 


Property changes and shrinkage effects at the laminate 
level are calculated by using the modified Classical Lami- 
nated Plate Theory (CLPT) (K.J. Bowles, NASA Lewis Re- 
search Center, personal correspondence). Laminate properties 
are calculated according to the relations in appendix B, with 
the altered ply properties. The effects of shrinkage on the 
laminate are analogous to those caused by thermal deforma- 
tions. A modification of CLPT is used to define and apply 
effective in-plane and bending degradation forces to the lami- 
nate. Next, laminate deformations are calculated, and then ply 
stresses are calculated from the deformed laminate shape. 

The shrinkage strain for each ply i, with ply angle 0., are 
expressed in compact form as 


k = b r (N + N r + N // + N*) + d(M + M r + M H + M*j (71) 

The a, b, and d matrices are the inverse laminate stiffness 
matrices, calculated according to the relations given in appen- 
dix B. The midplane strains in each layer are 

e«0 =e o +Kz < (72) 

where z- is the distance from the midplane of ply i to the 
laminate midplane. The stresses in the laminate coordinate 
system are then 



The shrinkage strains in the laminate coordinate system are 





(67) 


CT /(i) - Q;( e /(i) + £ t(i) + e Hi) + e ((i ) ) 

where £.J is the ply free thermal strain in the laminate coordi- 
nate system and is the ply free hygral strain in the lami- 
nate coordinate system. In the ply coordinate system 

l 

°m ■ =T i°m (74) 

°m L 
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7.2 Ply Failure 


8.2 Material Data Base 


The maximum stress failure criterion ratios for longitudi- 
nal, transverse, and shear stress in each ply are calculated in 
this section: 


R,= 


u n\ 

s n\T 

s n\c 


a m >0 


°m <0 


(75) 


R. 


u >22 
S n2 T 

s mc 


a (22 ^ 0 


CJ £22 ^ ^ 


(76) 


R = 




(77) 


The existing ICAN data base was updated with properties 
for the C-6000 fiber and modified properties for the PMR-15 
matrix material. The fiber properties were taken from refer- 
ence 8; unavailable properties were scaled from the existing 
data base for the similar AS4 fibers. The existing PMR-15 
data base was modified as follows: (1) the shear strength was 
increased to reflect the ratio of shear strength to the ultimate 
tensile strength observed for virtually all resin materials; 
(2) the initial value of T was updated (ref. 2); and (3) degra- 
dation material constants were added. Cross-linking was as- 
sumed to cause no material property changes other than those 
due to an increase in T (which are extensive at high tempera- 
tures), and to cause no shrinkage or density change. Mass loss 
was, as a conservative baseline, assumed to cause only shrink- 
age (no void formation). This assumption contradicts some 
available evidence (ref. 7), so the percentage of mass loss due 
to void formation was later studied parametrically. The effec- 
tive properties of the material were assumed to decrease in 
linear proportion to mass loss. This is consistent with assum- 
ing that the remaining material has unchanged material prop- 
erties, which is supported by a small amount of evidence 
(ref. 15). The baseline data used are found in appendix D. 


If R . > 1 .0, the ply has failed in mode L 

The ultimate tensile and compressive strengths ( T and C 
subscripts) in the fiber direction and transverse direction ( 1 1 
and 22 subscripts) and the ultimate shear strength S^ l2s are all 
calculated at the moisture, temperature, and degradation con- 
ditions of the ply. This is done by using existing microstruc- 
tural relations (ref. 12) and matrix constitutive properties as 
modified by equations (4) and (6). 

8.0 Implementation 

8.1 Computer Code 

The method just outlined for calculating laminate deforma- 
tions, ply stresses, and microstresses is implemented as a 
module of the ICAN computer code. A complete description 
of this code is given in reference 11. The user specifies 
laminate geometry, loadings, temperatures, moisture contents, 
and degradation states. A typical input file is shown in appen- 
dix C. The user can also modify a material data base file 
containing constituent material properties and the constants 
relating degradation to material property changes. An ex- 
ample data base, containing properties for both the C-6000 
fiber and the PMR-15 matrix, is given in appendix D. The 
code output is illustrated in appendix E. A useful spinoff of 
this work has been the creation of an easy-to-use Macintosh 
version of the ICAN code. 


9.0 Verification and Example Uses 

The code was exercised to verify that its predictions could 
correlate with experimentally observed results. Two verifica- 
tion correlations were performed: (1) Predictions of changes 
in the high-temperature properties of PMR matrix composites 
were made on the basis of assumed cross-linking and then 
were correlated with available data. (2) Predictions of gross 
sample shrinkage were made by using the data on mass loss 
and growth of a surface degradation zone with the mass loss 
degradation model, and the predictions were correlated with 
available data. The verified code was then used to reach a 
better understanding of (1) the causes of cracking in unloaded 
specimens, (2) the response of aged specimens undergoing 
longitudinal bending tests, and (3) the response of aged speci- 
mens undergoing transverse bending tests. 

9.1 Cross-Linking Correlation 

The bending and shear strengths of a unidirectional sample 
were predicted at various temperatures as functions of the 
cross-linking metric a . The results were correlated with the 
data of Bowles (ref. 9) and Vannucci and Cifani (ref. 2). 
Room temperature results from both authors were used, but 
only the after-aging results from Bowles were used. He aged 
his specimens in nitrogen; Vannucci and Cifani aged theirs in 
air, so their results are complicated by surface degradation 
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effects. The parameter T ,in equation (1), defined as the glass 
transition temperature of a fully cross-linked matrix, was not 
available in the literature. In fact, it was noted in reference 9 
that severely aged specimens display no measurable glass 
transition up to the pyrolysis temperature of the matrix (above 
400 °C). Therefore, the parameter 7^ was determined by fit- 
ting equations (1) and (6) to data for unaged (a = 0 assumed) 
and severely aged (a ~ 1 assumed) specimens. This resulted in 
T , = 620 °C, which is not physically meaningful, being well 
above the pyrolysis temperature of the matrix, but it does 
provide a reasonable fit to available data, as illustrated in 
figures 4 and 5. 

9.2 Mass Loss and Shrinkage Correlation 

Bowles et al. (ref. 7) provided complete information on the 
mass loss and shrinkage of neat resin specimens undergoing 
isothermal oxidative degradation at a number of different 
temperatures. They noted growth of a degradation layer near 
the surface and plotted thickness of this layer as a function of 
exposure time (ref. 7). These data are sufficiently detailed to 
validate the shrinkage model developed here. 

In this report, the mass-loss data and the measured thick- 
nesses of degraded layers are used to define the degradation 
state of the specimen in terms of the metrics defined in sec- 
tion 5. From the relations in sections 5 and 6, the response of 
the degraded layer is found. Next, the relations in section 7 
are used to predict the shrinkage of the entire specimen. The 
predicted and measured specimen shrinkages are then com- 



0.00 0.25 0.50 0.75 1.00 

Degree of matrix cross-linking 

Figure 4. — Longitudinal bending strength correlation. Data of 
Bowles (ref. 9) and Vannucci and Cifani (ref. 2). 


Prediction Experimental Temperature, 
data °C 

150 — O 23 

O 260 



0.00 0.25 0.50 0.75 1.00 

Degree of matrix cross-linking 

Figure 5. — Short beam shear strength correlation. Data of 
Bowles (ref. 9) and Vannucci and Cifani (ref. 2). 

pared. As an additional check on the method, the calculations 
(which are tractable due to the simple geometry of the prob- 
lem) were also carried out directly and checked against the 
code results. 

The specimens were 6.4 mm square and 75 mm long. For 
the purpose of this correlation, they are modeled as three- 
layer laminates — two degraded layers and a nondegraded cen- 
tral core (see figs. 6 and 7). Such a model is accurate for axial 
stress and shrinkage predictions, which depend solely on the 
relative areas of the degraded and undegraded material. Here, 
all mass loss was assumed to take place in the degraded 



Figure 6. — Cross section of aged square cross-section 
specimen showing degraded layer. 
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Figure 7. — Model of specimen as plate with same proportion 
of degraded material. 


layers, and within these layers the mass loss was assumed 
uniform. Therefore, in the core b - 0, and in the degraded 
layers 



(78) 



Exposure time, hr 


where Am is the measured mass change of the specimen and 
m d is the original mass of the material in the degraded layer as 
calculated from the geometry of the specimen, the density of 
the undegraded material, and the measured thickness of the 
degraded layer. 

The measured degradation thickness and the calculated 
degradation parameter were used in the model shown in fig- 
ure 7 to predict the longitudinal shrinkage of the specimen. 
Because of the changing thickness of the degraded layer, a 
new model was constructed for each temperature and expo- 
sure time. The shrinkage of the specimen was predicated on 
the assumptions that there were no mechanical loads or resid- 
ual thermal stresses and that a = 1 (based on the long exposure 
times), D = 0.1 (based on the observation of moderate void 
formation in the degraded layer), and the room temperature 
material properties were unchanged by mass losses. Paramet- 
ric studies showed that the results were relatively insensitive 
to these assumptions. Predicted and measured results are 
compared in figure 8. As can be seen, the correlation is quite 
good. 

The predicted results also provide insight into the response 
of these specimens during exposure. The mass loss in the 
surface layer results in shrinkage, which is constrained by the 
unshrinking, undegraded core material. The results are a net 
shrinkage of the entire specimen, and stresses occurring in 
both the degraded layers and the core. The stresses in the core 
are compressive, and the stresses in the degraded layer are 
tensile. These stresses reach a magnitude that is sufficient to 
cause failure in the degraded layer. Although insufficient data 


Figure 8. — Comparison of predicted and measured changes 
in longitudinal dimension (predictions are based on de- 
graded layer thickness and mass-loss data). 


were provided in reference 7 for a quantitative correlation of 
this effect, the predicted failure in the surface layers qualita- 
tively agrees with the observed cracking of the surface layer. 

93 Surface Cracking Study 

The consequences of mass loss near the surface of compos- 
ite specimens was parametrically studied. Figure 9 shows a 
[0/±45/90]s laminate with mass-loss degradation occurring 
only in the surface plies. The specimen was assumed to have 
been exposed to high-temperature conditions long enough to 
fully cross-link the matrix (a = 1). The stresses in the surface 
plies of an unloaded laminate were examined as mass-loss 
parameter b was increased. 

Figure 10 shows the fiber-direction and transverse stresses 
in the outer plies of the laminate shown in figure 9 as func- 
tions of mass loss in these plies. A constant temperature of 
371 °C is assumed. Large stresses occur in these plies because 
of the matrix shrinkage associated with mass loss. A trans- 
verse ply tensile strength of about 30 MPa is predicted at this 
temperature, implying the stresses are sufficient to cause ply 
cracking at about 5-percent mass loss in the surface plies. The 
maximum matrix stresses at the micromechanical level. 
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Matrix stress, MPa Ply stress, MPa 


300 . — 


Degraded 
outer plies 


[0/90/±45] s layup 


Figure 9. — Laminate with degraded surface plies. 



1501 - 



Mass loss, percent 

Figure 10. — Ply stresses in surface plies at 371 °C. 



Mass loss, percent 

Figure 11 . — Shrinkage-induced microstresses in matrix of 
surface plies at 371 °C. 



0 5 10 15 20 

Mass loss, percent 

Figure 1 2. — Ply stresses in surface plies at 25 °C. 


calculated from equations (28) and (30), are shown in fig- 
ure 1 1 . They are somewhat higher than the ply stresses. The 
ply stresses at room temperature are shown in figure 12; they 
are much higher than the stresses at high temperature because 
of the presence of residual curing stresses (the cure tempera- 
ture was assumed to be 371 °C) and the fact that the matrix is 
stiffer at room temperature. Although the matrix is also much 
stronger at room temperature, the stresses are easily sufficient 
to cause cracking in this case. Even at zero mass loss, the 
stresses are very close to the transverse strength of the matrix. 
This raises the possibility that specimens with small mass 
losses may not crack at the aging temperature, but may crack 
when returned to room temperature. 

Figures 13 to 15 show how varying other parameters effects 
the outer ply stresses. A baseline case of 371 °C and 5-percent 
mass loss in the outer plies was used. Figure 13 shows that if 
some mass loss is occurring in the bulk of the material (the 
default case assumed all of the mass loss was in the surface 
plies) the fiber-direction stresses are lowered, but not the 
critical transverse stresses. Figure 14 shows that if some of 
the mass loss is due to void formation instead of shrinkage, 
stresses are reduced proportionally. On the other hand, differ- 
ent assumptions about the reduction of room temperature 
material properties with degradation have little effect on the 
calculated stresses. Figure 15 shows the effect of different 
material degradation parameters A E t A 5 , and A a a (trans- 
verse properties), which varied between 0 and -1 .0, and their 
respective degradation exponents n £a , n Sa , and which 
varied from 0.5 to 1.0. 
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Ply stress, MPa 31 Ply stress, MPa 



gure 13. — Effect of bulk mass loss on surface ply stresses 
(temperature, 371 °C; surface ply mass loss, 5 percent). 



Property loss parameter, percent 


Figure 1 5. — Effect of material property changes on surface 
ply stresses (temperature, 371 °C; surface ply mass loss, 
5 percent). 



Figure 14. — Effect of void formation on surface ply stresses 
(temperature, 371 °C; surface ply mass loss, 5 percent). 


9.4 Longitudinal Bending Test Simulation 

A 16-ply unidirectional specimen 2.24 mm thick, aged at 
371 °C, was examined. Aging and degradation of the speci- 
men were assumed to be functions of time. The cross-linking 
metric a in all plies and the degradation metric b in the surface 
plies (b = 0 elsewhere) are shown in figure 16. These curves 
are approximations based on experimentally observed behav- 
ior. The assumption of rapid completion of cross-linking was 
based on information provided in reference 9. The mass loss 
was computed from equation (78) by using the assumed curve 
for b\ the results compare well with the data of reference 2. 
Note that the massive mass loss in the surface plies correlates 
with the observation that bare fibers eventually flaked off the 
specimens. 

The bending strength of the modeled specimens was com- 
puted by applying a trial bending load and comparing the 
stress in the outer ply on the compressive side of the specimen 
to the ultimate bending strength for this ply as computed by 
the code. The calculated results were iterated until the bend- 
ing load necessary to just fail the specimen was found. The 
results are presented in figure 17. Two types of “failure stress” 
are shown: the stress in the center of the outer ply at failure as 
calculated by the code, and the equivalent stress calculated for 
a uniform piece of material being tested. The data from refer- 
ence 2 were assumed to have been reduced in this way. This 
equivalent stress is 
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Maximum stress at failure, MPa Degradation metric 



Figure 1 6. — Assumed degradations as function of aging 
time. 



Figure 17. — Longitudinal bending test correlation. 



Stress, MPa 

Figure 18. — Comparison of ply stresses predicted by code 
to stresses calculated by assuming uniform material at 
time of failure for specimens aged 1 50 hr. 


a 
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where M is the applied bending load at failure, and H is the 
laminate thickness. The results are seen to be very good in 
terms of matching measured and calculated equivalent stresses, 
but the equivalent stress and the actual failure stress are not 
the same because of the nonuniformity of the material. This 
problem is illustrated in figure 1 8, which shows the stresses at 
the point of failure in each ply in a sample aged for 150 hr. 
The apparent stepwise variation in the stresses is an artifact of 
the code’s calculation of a single stress for each ply and is not 
important. However, the top and bottom plies have stresses 
that differ significantly from those predicted for a uniform 
material. This effect complicates accurate interpretation of 
bending test data. 


9.5 Transverse Bending Test Simulation 


The data of transverse bending tests have been found to be 
very sensitive to aging (ref. 16), which suggests that these 
tests may be useful for screening materials. A transverse 
bending test simulation was attempted with the code, but the 
test failed, for interesting reasons. It was found that under the 
same aging conditions assumed in the longitudinal bending 
simulation described in section 9.4, the surface plies cracked 
and lost transverse load-carrying capability, even without a 
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load being applied. If the remaining strength of the specimens 
were calculated by assuming the surface plies did not exist, 
the results were not conservative. Evidently the surface plies 
provide starter cracks that badly weaken the specimens. Such 
cracks have been noted, as has their tendency to become self- 
perpetuating since they provide a pathway for air and promote 
additional oxidation. This phenomenon is beyond current ana- 
lytical capabilities, and it greatly complicates the quantitative 
interpretation of the results of transverse bending tests. 

10.0 Conclusions 

A model of the effects of material degradation on the stress 
state in polymer matrix composites has been developed and 
demonstrated. It successfully replicates the behavior observed 
in laminates undergoing isothermal aging. The increase in 
high -temperature properties in the early phases of high- 
temperature exposure was assumed to be due to matrix cross- 
linking and was modeled semi-empirically. Further exposure 
effects were primarily mass loss and resultant shrinkage near 
the surface. Micromechanical and laminate theory relations 
for predicting ply and laminate behavior due to shrinkage and 
property changes were developed and verified. It was found 
that nonuniform shrinkage causes stresses, which can in turn 


cause ply cracking. Degradation also creates a nonuniform 
material, which makes the interpretation of existing test data 
difficult. The stress at failure in a longitudinal bending test is 
not the same as that calculated by assuming a uniform mate- 
rial, for example. The cracks also strongly affect the results of 
transverse bending tests, thereby making the results of these 
tests very difficult to use quantitatively. 
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11.0. Appendix A 
Symbols 


A laminate in-plane stiffness matrix 

a inverse laminate in-plane stiffness matrix 

a cross-linking metric 

B laminate coupling stiffness matrix 

b inverse laminate coupling stiffness matrix 

b mass-loss metric 

C shrinkage constant for cross-linking mechanism 

D laminate bending stiffness matrix 

D fraction of mass loss due to changes in porosity 

d inverse laminate bending stiffness matrix 

E Young’s modulus 

AE normalized change in matrix stiffness properties 

G shear modulus 

H laminate thickness 

h degraded layer thickness 

k volume fraction 

/ side length of control volume 

M applied moment vector 

M applied bending load at failure 

AM percent moisture content by weight 

m mass in control volume 

Am measured mass change 

N applied force vector 

N number of plies in laminate 

P(T) matrix property at temperature T 


Q reduced stiffness matrix 

Q rotated reduced stiffness matrix 

Rf failure stress ratio 

S composite strength 

AS normalized change in matrix strength properties 

T transform matrix 

T temperature 

T sf stress-free temperature of laminate 

AT change in temperature 

V volume of material in control volume 

z distance from centerline of ply to laminate centerline 

z? distance from bottom of ply to laminate centerline 

z l distance from top of ply to laminate centerline 

a thermal expansion coefficient 

Aa normalized change in matrix coefficient of thermal 

expansion 

moisture expansion coefficient 
8 specimen length change 

e strain vector 

i rotated strain vector 

£ 0 laminate strain vector 

e matrix strain 

£ shrinkage coefficient for cross-linking mechanism 

6 ply angle 

k laminate curvature vector 
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V 


Poisson’s ratio 


p 

density 

a 

rotated stress vector 

a 

stress in component 

cfq 

maximum stress in uniform beam at failure 

<t> 

porosity 

X 

shrinkage coefficient for mass-loss mechanism 

Superscripts: 

b 

bottom 

H 

hygral 

n 

shrinkage exponent for cross-linking mechanism 

n Ea' n Eb 

exponents of change in matrix stiffness properties 
associated with first and second degradation 
mechanisms 

n Sa' n Sb 

exponents of change in matrix strength properties 
associated with first and second degradation 
mechanisms 

n aa* n ab 

exponents of change in matrix coefficients of ther 
mal expansion associated with first and second 
degradation mechanisms 

T 

thermal 

T 

matrix transpose 

t 

top 

* 

shrinkage 

Subscripts: 

A 

region A in composite microstructure model 

a 

associated with first degradation (cross-link) 
mechanism 


B region B in composite microstructure model 

b associated with second degradation (mass- loss) 

mechanism 

C compressive 

d initially degraded 

/ in fiber direction 

g glass transition 

gf final glass transition 

i relating to ply / 

ij in ij direction 

€ composite 

€\ 1 composite fiber direction 

( 1 2 composite shear direction 

€12 composite transverse direction 

m matrix 

o initial value 

RT room temperature 

S shear 

s shear 

T tensile 

t in transverse direction 

v void 

x component x 
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13.0 Appendix B 

Summary of Classical Laminate Theory of Equations 


For each ply i, we have ply properties E em (fiber-direction 
stiffness), £^ 22 - ( transverse stiffness), v^ 12j . (major Poisson’s 
ratio), Gf l2i (shear stiffness), a €Ui and $ €ni (fiber direction 
coefficient of thermal expansion (CTE) and moisture swelling 
coefficient), and a^ 22/ and fi n2i (transverse CTE and coeffi- 
cient of mechanical expansion (CME)). The fibers are aligned 
at an angle to the structural jc-axis. For each ply i, 

@11 (/) Q\1 (i) 0 

Q/ = @\2(i) @22(i) ° 

0 0 266(i) 

where 

Qim =E mi/ D i 

G 12(i) = v (l 2 i E e 22 i/ D i 



2 


and 




D 4i>[(*;) 3 4?) 3 ] lB3) 

1=1 

where z \ and Z* represent the distances from the top and the 
bottom of the ply to the laminate centerline, respectively. 
Assembling these 3x3 matrixes into a 6x6 matrix and invert- 
ing yields 



e 22 (i) “ E 122 i/ D i 
Q 66 (i) = G e\2i/ D i 

D (i) = 1 + v ll 2i E i22i l E l\ 1 i 

and 

Qf = T^Q.Tr 1 (B2) 

for 

cos 2 0. sin^ 2 sin 0- cos 0. 

T. = sin 2 0 / cos 2 0. -2 sin 0. cos0. 

2 2 

-sin 0. cos0. sin 0. cos 0- cos 0. - sin 0- 
For a laminate with N plies. 


where a, b, and d are the compliance submatrices for the 
laminate. 

For AT. = T. - where 7^is the stress-free temperature of 

the laminate, T. is its temperature, and AM is its percent 
moisture content by weight, the ply thermal strains are given 
by 



The ply moisture strains are given by > (B5) 

'ftl 

4)=^'AA/,. 

0 . 
i 

and the corresponding rotated ply strains about the structural 
x-axis are given by 
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and the laminate moments are given by 
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13.0 Appendix C 
Modified Input File 


A complete description of the input to the ICAN code is 
contained in reference 11. Here we will illustrate the new 
input features by an example input file, with new input fields 
in bold. This file must reside in the same folder as ICAN. It 
can have any name; ICAN will ask for an input file name 
when it is launched. 


Parametric Study input file, C-6000/PMR-15 composite 
COMSAT T 
CSANB F 
BIDE F 
RINDV F 
NONUDF T 
DEFECT F 


PLY 

1 

2 

700.00 

700.00 

0.00 

0.000.005500 

1.00 

0.10 

PLY 

2 

1 

700.00 

700.00 

0.00 

45.000.005500 

1.00 

0.01 

PLY 

3 

1 

700.00 

700.00 

0.00 

-45.000.005500 

1.00 

0.01 

PLY 

4 

1 

700.00 

700.00 

0.00 

90.000.005500 

1.00 

0.01 

PLY 

5 

1 

700.00 

700.00 

0.00 

90.000.005500 

1.00 

0.01 

PLY 

6 

1 

700.00 

700.00 

0.00 

-45.000.005500 

1.00 

0.01 

PLY 

7 

1 

700.00 

700.00 

0.00 

45.000.005500 

1.00 

0.01 

PLY 

8 

2 

700.00 

700.00 

0.00 

0.000.005500 

1.00 

0.10 

MATCRD 

1C-6KPMHM 0.5800 

0.000 1C-6KPMHM 

0.0000 0.0000 

0.0000 


MATCRD 

2C-6KPMHM 0.5800 

0.0001 C-6KPMHM 

0.0000 0.0000 

0.0000 



PLOAD 1000. 
PLOAD 
PLOAD 
OPTION 5 
OPTION 7 
OPTION 1 1 
OPTION 12 
OPTION 15 


The two additional eight-character fields on the PLY cards 
are the input values of the degradation metrics a and b for 
each ply. In this example, all plies are fully cross-linked, the 
surface plies have suffered 10-percent matrix mass loss, and 
the inner plies have suffered 1 -percent matrix mass loss. The 
existing ICAN convention, which requires a unique MATCRD 
card for plies with different temperatures, applies to the deg- 
radation state. For each unique degradation state, a separate 
MATCRD card must be used. In this example, two states 
exist, distinguished by two different values of b , so two 
MATCRD cards are used. 

No choice of output OPTION is required. If any nonzero 
degradation is specified, the degradation-induced property 
and stress output will be provided automatically. 
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14.0 Appendix D 
Degraded File 


A complete description of the material files used by the 
ICAN code is contained in reference 1 1 . Here we list the new 
material properties used for this study and illustrate the new 
degradation constant inputs. Material properties are stored in 
a file named MAT.DAT, which must reside in the same folder 
as the ICAN code. Degradation constants for matrix materials 
are preceded by a single line containing only the characters 
DEGRADATION CONSTANTS. This line must immedi- 
ately follow the regular material properties. If degradation 
constants are not provided, but nonzero degradation is speci- 
fied for any plies, then default values are used. 


FIBER PROPERTIES 

C-6K C-6000 EQUIVALENT FIBER PROPERTIES (H. MCMANUS 8-93) 
$ 

$ 

$ 


Number of fibers per end 

Nf 

15000 

number 

Filament equivalent diameter 

df 

0.280E-03 

inches 

Weight density 

Rhof 

0.640E-01 

lb/in* *3 

Normal moduli (11) 

Efll 

0.340E+08 

psi 

Normal moduli (22) 

Ef22 

0.220E+07 

psi 

Poisson’s ratio (12) 

Nufl2 

0.220E+00 

non-dim 

Poisson’s ratio (23) 

Nuf23 

0.250E+00 

non-dim 

Shear moduli ( 1 2) 

Gfl2 

0.540E+07 

psi 

Shear moduli (23) 

Gf23 

0.110E4O7 

psi 

Thermal expansion coef. (11) 

Alfafl 1 

-590E-06 

in/in/F 

Thermal expansion coef. (22) 

Alfaf22 

0.560E-05 

in/in/F 

Heat conductivity (11) 

Kfll 

0.440E+01 

BTU-in/hr/in**2/F 

Heat conductivity (22) 

Kf22 

0.440E+00 

BTU-in/hr/in* *2/F 

Heat capacity 

Cf 

0. 1 66E+00 

BTU/lb/F 

Fiber tensile strength 

SfT 

0.550E+06 

psi 

Fiber compressive strength 

SfC 

0.650E+06 

psi 


OVER END OF FIBER PROPERTIES 


MATRIX PROPERTIES 

PMHM PMR-15 POLYIMIDE MATRIX modified by H. McManus 8-93 
$ increased initial Tg to match Vannucci data 
$ increased shear strength 
$ 


Weight density 

Rhom 

0.440E-01 

Ib/in**3 

Normal modulus 

Em 

0.470E+06 

psi 

Poisson’s ratio 

Num 

0.360E+00 

non-dim 

Thermal expansion coef. 

Alfa m 

0.280E-04 

in/in/F 

Matrix heat conductivity 

Km 

0.146E+00 

BTU-in/hr/in**2/F 

Heat capacity 

Cm 

0.250E+00 

BTU/lb/F 

Matrix tensile strength 

SmT 

0.810E+04 

psi 

Matrix compressive strength 

SmC 

0.165E+05 

psi 

Matrix shear strength 

SmS 

0.135E+05 

psi 

Allowable tensile strain 

eps mT 

0.183E-01 

in/in 

Allowable compr. strain 

eps mC 

0.350E-01 

in/in 

Allowable shear strain 

eps mS 

0.500E-01 

in/in 

Allowable torsional strain 

eps mTOR 

0.500E-01 

in/in 

Void heat conductivity 

kv 

0.225E+00 

BTU-in/hr/in* *2/F 

Glass transition temperature 

Tgdr 

0.740E+03 

F 

DEGRADATION CONSTANTS 



Max aged glass transition T 

Tgmax 

0.1 15E+04 

F 

Shrinkage due to crosslinking 

Cal 

O.OOOE+OO 

non-dim 

exponent for above 

nal 

1.000E+00 

non-dim 

void frac delta from mass loss 

D 

0.000E+00 

non-dim 

Changes in properties (D) 

DEI 

0.000E+00 

non-dim 

(dimensionless fraction of 

nEl 

1.000E+00 

non-dim 

original properties) 

DE2 

100E+01 

non-dim 

and corresponding exponents 

nE2 

1.000E+00 

non-dim 

(n) for crosslink (1) and 

DS1 

0.000E+00 

non-dim 

massloss (2) aging 

nSl 

1.000E+00 

non -dim 

E = modulus 

DS2 

100E+01 

non -dim 

S = all strengths 

nS2 

1.000E+00 

non-dim 

alf = CTE 

Dalfl 

O.OOOE+OO 

non -dim 


nalfl 

1.000E+00 

non-dim 


Dalf2 

1 00E+0 1 

non-dim 


nalf2 

1.000E+00 

non-dim 
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15.0 Appendix E 
Output File 


The following pages contain an abbreviated version of an 
output from the ICAN code run by using the input file in 
appendix C and the material file in appendix D. Omitted 
material, which is unchanged from the original ICAN, is 
indicated by three stars, i.e., ***. 

input file name? 

ninput «< SERE HE PROVIDE THE NAME OF THE INPOT FILE 


Editorial comments added to the new output are in bold be 
tween angles 

<«LIKE THIS>» 

INTERACTIVELY »> 


ICAN 

INPUT DATA ECHO 

* ** 


SUMMARY OF INPUT DATA 
Parametric Study input file, C6000/PMR15 composite 

** * 

- * - LAMINATE CONFIGURATION - - - <« NOTE DEGRADATION METRICS »> 


PLY NO MID DELTAT DELTAM THETA T-NESS C-LINK D-GRADE 


PLY 

1 

2 

0.000 

0.0% 

PLY 

2 

1 

0.000 

0.0% 

PLY 

3 

1 

0.000 

0.0% 

PLY 

4 

1 

0.000 

0.0% 

PLY 

5 

1 

0.000 

0.0% 

PLY 

6 

1 

0.000 

0.0% 

PLY 

7 

1 

0.000 

0.0% 

PLY 

6 

2 

0.000 

0.0% 


0.0 

0.005 

1.000 

0.100 

45.0 

0.005 

1.000 

0.010 

-45.0 

0.005 

1.000 

0.010 

90.0 

0.005 

1.000 

0.010 

90.0 

0.005 

1.000 

0.010 

-45.0 

0.005 

1.000 

0.010 

45.0 

0.005 

1.000 

0.010 

0.0 

0.005 

1.000 

0.100 


*** 


--> CONSTITUENT PROPERTIES: ECHO FROM DATA BANK . 
PRIMARY FIBER PROPERTIES; C-6K FIBER 


1 

ELASTIC MODULI 

EFP1 

0.3400E+08 

2 


EFP2 

0 . 2200E+ 07 

3 

SHEAR MODULI 

GFP12 

0 . 54 00E + 07 

4 


GFP23 

0 . 1 100E+07 

5 

POI SSON M S RATIO 

NUFP12 

0 .2200E+00 

6 


NUFP23 

0 .2500E+00 

7 

THERM. EXP. COEF. 

CTEFP1 

-0 . 5900E- 06 

8 


CTEFP2 

0 . 5600E-Q5 

9 

DENSITY 

RHOFP 

0 . 6400E-01 

10 

NO. OF FIBERS/END 

NFP 

0 . lQOQE-t-0 1 

11 

FIBER DIAMETER 

DIFP 

0.2800E-03 

12 

HEAT CAPACITY 

CFPC 

0 , 1660E+00 

13 

HEAT CONDUCTIVITY 

KFPl 

0 .4400E+01 

14 


KFP2 

0 .4400E+00 

15 


KFP3 

0 .4400E+00 

16 

STRENGTHS 

SFPT 

0 . 5500E+06 

17 


SFPC 

0 . 6500E+06 


«< HERE HE COMPUTE THE DEGRADED PROPERTIES. NOTE ECHO OF DEGRADATION CONSTANTS >» 


PRIMARY MATRIX PROPERTIES ; PMHM MATRIX. DRY RT . PROPERTIES. 


REFERENCE TEMP. - 70.00 

CROSSLINK METRIC - 1.000 

ORIG. DRY GLASS TRANS. TEMP. 
DEGRADATION CONSTANTS FOR THIS 
TGMAX - U50.000 CA1 - 

DEI - 0.000 NE1 - 

DS1 - 0.000 NS1 - 

DALF1 - 0.000 NALF1 - 


TEST 

TEMP. 

700.00 

PCT . 

MOISTURE 

DEGRADATION 

METRIC - 0 . 

010 


740.00 


MODIFIED BY 

CROSSLINKING 

- 11 SO. 

MATERIAL - 
0.000 

NA1 

1.00C 

CA2 - 

0.000 

1.000 

DE2 

-1.000 

NE2 - 

1.000 

1 .000 

DS2 

-1.000 

NS2 - 

1.000 

1,000 

DALF2 

-1.000 

NALF2 - 

1.000 


0.000 


ORIGINAL MATRIX PROPERTIES 


1 ELASTIC MODULUS 

2 SHEAR MODULUS 

3 POISSON ,f S RATIO 

4 THERM. EXP. COEF. 

5 DENSITY 

6 HEAT CAPACITY 

7 HEAT CONDUCTIVITY 


EMP 

D.4700E+06 

GMP 

0. 1728E+06 

NUMP 

0 . 3600E+00 

CTEMP 

0.2800E-04 

RHGMP 

0.4400E-01 

CMPC 

0 . 2500E+00 

KMP 

0. 1460E + 00 


NEW PROPERTIES: CHANGED BECAUSE 
CF TEMP /MOIST/CLINK/DGRD EFFECTS 


NEW EMP 
NEW GMP 
NEW NUMP 
NEW CTEMP 
NEW RHOMP 
NEW CMPC 
NEW KMP 


0.3003E+06 
0. H04E+06 
0.3600E+00 
0 . 4294E-04 
0 .4400E-01 
0 . 3873E+00 
0 . 2262E+ 00 
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STRENGTHS 


9 

10 

11 MOISTURE COEF 

12 DIFFUS IVITY 


SMPT 0-81 00E+04 
SMPC 0.16S0E+05 
SMPS 0 - 1 350E + 0 5 
BTAMP 0 . 4 OOOE-02 
DIFMP 0 . 2D00E- 03 


NEW SMPT 
NEW SMPC 
NEW SMPS 
NEW BTAMP 
NEW DIFMP 


O.S176E+04 
0 . 10S4E+0 5 
0 .8627E + 04 
0.4000E-02 
0.20Q0E-03 


*** «< A LARGE AMOUNT OF OUTPUT OMITTED HERE- DEGRADED PROPERTIES FOR MATERIAL 2, COMPOSITE PROPERTIES, 
AND THE STRESS CALCULATIONS FROM THE ORIGINAL I CAN (WHICH INCLUDE DEGRADATION- INDOC ED CHANGES 
IN THE MATERIAL PROPERTIES, BUT NOT THE EFFECTS OF MATERIAL SHRINKAGE) »> 


«< NOW WE OUTPUT LAMINATE DISPLACEMENTS AND PLY STRESSES DUE SOLEY TO DEGRADATION >» 


LAMINATE DISPLACEMENT-FORCE RELATIONS SOLELY DUE TO EFFECTS OF MATRIX DEGRADATION 


DISP. 


DISPIACEMENT FORCE RELATIONS 


DEGRAD. FORCES 


1 

2 

3 

4 

5 

6 


-0. 1036E-O3 
-0 . 2906E-03 
-0 . 2485E-1 0 
0 . 2594E-09 
0.9368E-10 
0. 7861E-09 


- 1 - 

0.3131E-05 
-0.9951E-06 
-0 . 8044E-13 
0.1007E-10 
-0.U53E-10 
-0.3427E-11 


- 2 - 

-0 . 993 IE- 06 
0.3135E-05 
0 . 2 S34E-12 
-0.4668E-11 
0.1066E-10 
-0.6329E-11 


-3- 

-0 . 8Q44E-I3 
0 .2534E-12 
0 .8277E-05 
-0.1819E-12 
-0.2I27E-11 
-0 .3200E-10 


-4- 

0.1007E-10 
-0.4668E-11 
-0.1819E-12 
0.1138E-Q1 
-0.7273E-02 
-0 . 1596E-02 


-5- 

-0. 1153E-10 
0. 1D66E-10 
-0.2127E-11 
-0.7273E-02 
0.5534E-01 
-0. 1866E-01 


- 6 - 

-0.3427E-11 

-0.6329E-11 

-0.3200E-10 

-0.1596E-02 

-0.1866E-01 

0.6915E-01 


-0. 9797E+02 
-0. 1238E+03 
-0. 1642E-06 
0.6786E-07 
0 . 1423E-07 
0.S898E-09 


NOTE: THE DISPLACEMENTS ARE REFERENCE PLANE MEMBRANE STRAINS <UX , VY , VXPUY) AND CURVATURES {WXX , WYY , WXY ) 


PLY STRESSES, STRAINS, AND MAX. STRESS CRITERIA FAILURE RATIOS 
DUE SOLELY TO DEGRADATION EFFECTS 
IN PLY {MATERIAL) COORDINATE SYSTEM 


PLY 

SIGH 

SIG22 

SIG12 

EPS11 

EPS22 

EPS 12 

SRI L 

SR 22 

SR 12 

1 

2 . 590E+03 

8.643E+03 

-3 . 861E-06 

-1 .836E-04 

-2. 906E-04 

-9.717E-12 

0.0081 

2.2227 

0.0000 

2 

-4 . 093E+03 

6 . 882E+02 

-4 . 650E+01 

*2 . 37 IE-04 

-2.371E-04 

-1.070E-04 

0.0497 

0.1605 

0.0067 

3 

-4 . 093E+03 

6 . 882E+02 

4 . 650E + 01 

-2.371E-04 

-2-371E-04 

1.07GE-04 

0.0497 

0.1605 

0.0067 

4 

-5 . 145E+D3 

7.221E+02 

1 .393E-05 

-2.906E-04 

-1 .836E-04 

3 . 204E-11 

0.0624 

0. 1684 

0.0QD0 

5 

-S.145E+03 

7 .221E+02 

1 .581E-05 

~ 2 . 906E-O4 

-1.836E-04 

3.636E-11 

0.0624 

0.1684 

C . 0000 

6 

-4.093E+03 

6.882E+02 

4 . 650B + 01 

-2 . 3 7 IE- 04 

-2.371E-04 

1 .070E-04 

0.0497 

0. 1605 

C . 0067 

7 

-4 . G93E + Q3 

6.882E+02 

-4 . 650E+Q1 

-2 .371E-04 

-2.371E-04 

-1 .070E-D4 

0.0497 

0.1605 

0.0067 

8 

2 . 590E+Q3 

8 .643E+03 

*1 .589F.-05 

•1 .836E-04 

-2. 906E-04 

-3.998E-11 

0.0081 

2.2227 

C . 0000 

«< NOW (FINALLY) WE OUTPUT 

THE DISPLACEMENTS AND 

STRESSES DUE TO THE ENTIRE LOADING 

CASE SPECIFIED- 


MECHANICAL, THERMAL, HYGRAL, AND DEGRADATION INDUCED »> 


LAMINATE DISP LACE. ME NT “FORCE RELATIONS INCLUDING ALL LOADING EFFECTS 
DISP. DISPIACEMENT FORCE RELATIONS COMBINED FORCES 




-1- 

-2- 

“3- 

-4- 

-5- 

-6- 


1 

0.2948E-02 

0.3131E-05 

0.9951E-06 

-0.8044E-13 

0. 1007E-10 

-0 . 1153E-10 

-0 . 3427E- 1 1 

0 . 9D20E+03 

2 

-0 . 1286E-02 

-0 . 9951E-06 

0. 3135E-05 

0 . 2534E-12 

-0.4668E-11 

0.1066E-10 

-0.6329E-11 

-0. 1238E+03 

3 

-0.10S3E-09 

-0.8044E-13 

0.2534E-12 

0.8277E-05 

-0.1819E-12 

-0.2127E-11 

-0.320QE-10 

-0. 1642E-06 

4 

0. 1033E-07 

0.1007E-10 

-0.4668E-11 

-0.1819E-12 

0. 1138E-01 

-0 . 7273E-02 

-0.1 596E-02 

0.6786E-07 

5 

-0. 1144E-07 

-0 . 1 153E-1 0 

0.1066E-10 

-0 . 2127E-11 

-0.7273E-02 

0 . 5534E-01 

-0.1866E-01 

0.1423E-07 

6 

-0 . 2641E-08 

-0.3427E-11 

-0.6329E-11 

-0.3200E-10 

-0 . 1 596E-02 

-0.1866E-01 

0 . 6915E-01 

Q.5898E-09 


NOTE: THE DISPLACEMENTS ARE REFERENCE PLANE MEMBRANE STRAINS <UX , VY , VXPUY) AND CURVATURES (WXX , WYY , WXY) 


PLY STRESSES, STRAINS, AND MAX. STRESS CRITERIA FAILURE RATIOS 
INCLUDING EFFECTS OF MECHANICAL LOAD, TEMPERATURE, MOISTURE, AND DEGRADATION 
IN PLY (MATERIAL) COORDINATE SYSTEM 


PLY 

SIGH 

SIG22 

SIG12 

EPS 1 1 

EPS22 

EPS 12 

SRI 1 

SR22 

SR12 

1 

6.467E+04 

8 . 542E+03 

- 6. 204 E- OS 

2 .948E-03 

-1.286E-03 

-1 .561E-10 

0.2027 

2.1968 

0.0000 

2 

1 . 744E+04 

1 .690E+03 

-1 . 840E + 03 

8 .311E-04 

6 .31 IE- 04 

4 . 233E-03 

0.0547 

0.4408 

0.2640 

3 

1 .744E+04 

1 .890E+03 

1 .840E+03 

8.311E-04 

8.311E-04 

4 . 233E-Q3 

0.0547 

0.4408 

0.2640 

4 

-2.419E+04 

3.234E+D3 

2.098E-04 

*1 .286E-03 

2.948E-03 

4.826E-10 

0.2936 

0.7541 

0.0000 

5 

-2.419E+04 

3 . 234E+03 

2.035E-04 

-1 .286E-03 

2. 948E-03 

4.681E-10 

0.2936 

0.7541 

0.0000 

6 

1 .744E+04 

1 . 690E+03 

1 . 840E + 03 

8.311E-04 

8.31 IE -04 

4 .233E-03 

0.0S47 

0.4408 

G . 2640 

7 

1 . 744E+04 

1 .890E+03 

-1.840E+03 

8.311E-04 

8.311E-04 

-4.233E-03 

0.0547 

0.4408 

0 . 2640 

8 

6.467E+04 

8 . 542E+03 

-2.163E-05 

2.948E-03 

-1.286E-03 

5 . 444E-11 

0.2027 

2.196B 

0.0000 


«< NOTE TRANSVERSE STRESS RATIO IN THE OUTER PLIES (SR22) IS > 1, INDICATING CRACKING »> 

«< FINALLY, COMPUTE MICROSTRESSES, FIRST DUE TO MICRO-LEVEL DEGRADATION EFFECTS ALONE, AND THEN 
COMBINED WITH THE EFFECTS OF DEGRADATION INDUCED PLY STRESSES. TO GET TOTAL MICROSTRESSES, 

YOU MUST COMBINE THESE WITH THE OTHER MICROSTRESSES CALCULALATED BY THE OLDER ROUTINES IN ICAN »> 


23 



MICROSTRESSES IN PRIMARY AND SECONDARY COMPOSITES 
DUE SOLELY TO DEGRADATION EFFECTS (NOT INCLUDING PLY STRESSES) 

PRIMARY VALUES ON FIRST LINE, SECONDARY (IF ANY) FOLLOWS 


PLY SMI1 


SF1 1 SM22A 5H22B SF22 SM33A SM33B SF33 SINCN 

+ + — — * + + + + + 


1 

9.049E+03 

-6 . 55 1E+03 

6. 1 54E+03 

-1 .927E+03 

-1 

2 

9. 94 8E+02 

“7 . 202E+02 

6 , 770E+Q2 

-2 . 120E+02 

-2 

3 

9. 948E+02 

-7.202E+02 

6 . 770E+02 

-2 . 1 20E+02 

-2 

4 

9. 948E+Q2 

-7 . 202E+02 

6 . 770E+Q2 

- 2 . 120E+02 

-2 

5 

9« 94 8E+02 

-7.202E+02 

6. 770E+02 

-2. 120E+02 

-2 

6 

9 . 948E+02 

-7.202E+02 

6.770E+02 

-2 . 120E+02 

-2 

7 

9. 94 8E+02 

-7.202E+02 

6.770E+02 

-2. I20E+02 

-2 

8 

9 . 049E+03 

-6 . 551E+03 

6 . 1 54E+03 

-1 .927E+03 

-1 


927E+03 

6. 1S4E+03 

-1 .927E+03 

-1 . 927E+03 

S.977E+03 

120E+02 

6.770E+02 

-2 . 120E+02 

-2 . 120E+O2 

6 . S76E+02 

120E+02 

6 . 770E+02 

-2.120E+02 

-2 . 120E+02 

6 . 576E+02 

I20E+02 

6. 77QE+02 

-2 . 120E+02 

-2. I20E+02 

6.576E+D2 

120E+02 

6.770E+02 

-2 . 120E+02 

-2.120E+02 

6.576E+02 

120E+02 

6. 770E+02 

-2.120E+02 

-2.120E+02 

6.576E+02 

120E+02 

6 . 77QE+02 

-2 . 120E+02 

-2.120E+02 

6.576E+02 

927E+03 

6. 154E+03 

-1.927E+03 

- 1 . 927E+03 

5.977E+03 


MICROSTRESSES IN PRIMARY AND SECONDARY COMPOSITES 
DUE TO DEGRADATION EFFECTS, INCLUDING PLY STRESSES 

PRIMARY VALUES ON FIRST LINE, SECONDARY (IF ANY) FOLLOWS 


PLY 

SMI 1 

SF1 1 

SM22A 

SM22B 

1 

9. 126E+03 

-2.072E+03 

9 . 579E+03 

8.349E+03 

2 

9.365E+02 

-7 . 730E+03 

9 . 515E+-02 

6.057E+02 

3 

9 . 365E+02 

-7 . 730E+03 

9 . 51SE+02 

6.057E+02 

4 

9.207E+02 

-9.533E+03 

9 . 640E+02 

6 . 464E+02 

S 

9. 207E+02 

-9 . 533E+03 

9 . 640E+02 

6 . 464E+02 

6 

9. 365B+02 

-7.730E+03 

9. S15E+02 

6 . 057E+02 

7 

9.365E+02 

-7 , 730E+03 

9. 515E+Q2 

6 . 057E+02 

8 

9. 126E+03 

-2 . 072E+03 

9 . 579E+03 

8.349E+03 


SF22 

SM33A 

SM33B 

SF33 

SINCN 

. 349E+03 

5 . 920E+03 

-1 . 854E+03 

-1 . 854E+03 

1 . 5S5E+04 

. 057E+02 

6. S26E+C2 

-2 . 044E+02 

-2 . 044E+C2 

1 . 418E+03 

. 057E+02 

6.526E+02 

-2.044E+02 

- 2 . 044E+02 

1 .418E+03 

. 464E+02 

6 . S04E+O2 

-2.037E+02 

-2.037E+02 

1.454E+03 

. 464E+02 

6. S04E+02 

-2.037E+02 

-2.037E+02 

1 . 454E+03 

. 057E+02 

6.526E+02 

-2 .044E+02 

-2 . 044E+02 

1.418E+03 

. 057E+02 

6. 526E+02 

-2.044E+02 

-2.044E+02 

1 . 418E+03 

. 349E+03 

5. 920E+03 

-1 . 854E+03 

-1 . 854E+03 

1 . 5S5E+04 


*** 

STOP 


SINCS 


O.OOOE+OO 
0 . OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 
O.OOOE+OO 


SINCS 


-9.801E+Q3 
-7 . 840E+02 
-7.840E+02 
-8 . 227E+02 
-8.227E+02 
—7 . 84 DE+02 
-7 . 84 QE+02 
-9 . 801E+03 
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